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ABSTRACT 



This paper continues a series in which we intend to show how all observables of galaxy clusters can be combined to recover the two- 
dimensional, projected gravitational potential of individual clusters. Our goal is to develop a non-parametric algorithm for joint cluster 
reconstruction taking all cluster observables into account. In this paper, we begin with the relation between the Compton-y parameter 
and the Newtonian gravitational potential, assuming hydrostatic equilibrium and a poly tropic stratification of the intracluster gas. 
We show how Richardson-Lucy deconvolution can be used to convert the intensity change of the CMB due to the thermal Sunyaev- 
Zel'dovich effect into an estimate for the two-dimensional gravitational potential. Synthetic data simulated with characteristics of the 
ALMA telescope show that the two-dimensional potential of a cluster with mass 5 x 10 14 h~ l M a at redshift 0.2 is possible with an 
error of < 5 % between the cluster centre and a radius r < 0.9 h' 1 Mpc. 
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1. Introduction 



Beginning with (Navarro et al. 1997 1996 hereafter NFW), a 
multitude of numerical simulations has shown that gravitation- 
ally bound structures dominated by dark matter should exhibit 
a universal density profile with three characteristic properties: 
it starts flat in the core, steepens around a scale radius r s , and 
asymptotically approaches a double-logarithmic slope near -3 
towards the virial radius r v ;, (e.g. Jing & Suto 2000; Merritt et al. 
2006 , |Navarro et al.|20041|Power et al.|20031|Moore et al.fT998T 



1999 see also Einasto & Haud 1989). The concentration param- 
eter c = r V vlr s is found in simulations to depend only weakly 
on the mass M. For cold dark matter, it decreases oc M~ 0A (see 



|Navarro et al.|1997|[T996{|Bullock et al.|2001||Dolag et al.|2004 



Duffy etal . 2008 1 Ekeetal.|2001||Gaoet al. 2008 ; M accio et al 



2008. |2007] |Neto et al.||2007[ |Seljak||2000t |Shaw et al.||2006[ 
Zhao et al.|2009 for examples) 



Moreover, cold dark matter is expected to clump on virtu- 
ally all scales. Massive objects such as galaxy clusters should 
thus have a broad spectrum of massive sublumps embedded (cf. 
Boylan-Kolchin et a l. 2009; Dolag et al. 2009; Gao et al.|20TT 



2004||Giocorietal.|2 010, Zentner et al. 2005) 



While these mostly qualitative statements hold for gravi- 
tationally bound structures composed of cold dark matter on 
all mass scales, the dark-matter distribution in galaxy clusters 
should be least affected by baryonic physics because in them, the 
cooling time of the baryonic matter almost everywhere exceeds 
the Hubble time. Galaxy clusters are thus perhaps the class of 
objects best suited for testing whether the expectations raised by 
simulations are indeed supported by real dark-matter structures. 



Galaxy clusters provide a multitude of observables. Weak 
and strong gravitational lensing effects allow direct inversions 
yielding the scaled surface-mass density or, equivalently, the 
effective lensing potential of a lensing mass distribution. X- 
ray emission and the thermal Sunyaev-Zel'dovich (hereafter 
SZ) effect reflect the physical state of the hot intracluster gas. 
Assuming equilibrium and stability, the gas properties can also 
be related to the gravitational potential. This suggests to devise 
a method by which lensing, X-ray, and thermal SZ data can be 
combined in a joint analysis aiming at recovering the gravita- 
tional potential best compatible with all observables. 

The advantage and at the same time the main obstacle in 
such an approach is that the observables probe the gravitational 
potential on different scales. Strong lensing typically occurs on 
angular scales smaller than » (20 - 30)", which at the same 
time approximately equals the resolution limit of weak lensing. 
Strong and weak lensing can be combined in the same analy- 
sis either by assuming a parametric form of the density pro- 
file to be adapted to both effects, or in a parameter-free ap- 
proach fitting the effective lensing potential on a grid whose 
resolution is adapted to the angular scales on which the effects 
occur. Based on a maximum-likelihood cluster-reconstruction 
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method (Bartelmann et al. 1996 1, we have in the past years devel- 
oped a multi-scale approach to parameter-free reconstructions of 
the lensing potential which has been shown to perform reliably 
with simulated data (|Cacciato et al.||2006[ |Merten et al.||2009[ 
Meneghetti et al. 2010), and which has been apphed to several 
clusters so far. The method adapts potential values in grid cells 
until the lensing data are reproduced optimally in a least-^ 2 sense 
(see also [Bradac et al.|2005l [20061 ) . 

The series of papers which the present paper belongs to aims 
at extending this method towards including all available cluster 
observables in a joint reconstruction aiming at the one gravita- 
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tional potential underlying all of them. In the first paper ( |Konrad| 
et al. 2013|l, we have shown how X-ray data can be analysed 



to recover the gravitational potential projected along the line- 
of-sight. Assuming hydrostatic equilibrium and a polytropic gas 
stratification, we could show that the projected gravitational po- 
tential can be reproduced with typical errors at the per-cent level 
under typical conditions. For simplicity, not for necessity, we 
have further assumed spherical symmetry there. Analysed this 
way, the X-ray data can now be combined with the lensing data 
by adding a suitable term to the x 2 function initally defined for 
lensing. 

This paper describes how this approach can be extended 
towards including thermal SZ data. Again for simplicity, we 
adopt a spherical galaxy-cluster model for simulating observa- 
tions with the signal and noise characteristics of ALMA. 

We begin in Sect. 2 with a brief overview of the thermal SZ 
effect and review the essential steps of our method. We keep this 
short and refer the reader to the detailed description in the first 
paper of this series where possible. In Sect. 3, we describe how 
we simulated thermal SZ data sets including a realistic noise 
model. We test our algorithm varying its input parameters and 
estimate the error of the reconstruction in Sect. 4. We estimate 
the quality of our reconstruction of the lensing potential by boot- 
strapping 200 realizations of one specific example of a galaxy 
cluster. In Sect. 5, we draw our conclusions. 



2. Recovering the projected gravitational potential 
from SZ data 

2. 1 . The thermal Sunyaev-Zel'dovich effect 



The thermal Sunyaev-Zel'dovich effect (SZ) ( Sunyaev &| 
|Zeldovich||1980| l is caused the hot electrons in the intracluster 
plasma that inversely Compton scatter the much less energetic 
photons of the cosmic microwave background (CMB) to higher 
energies. The SZ effect thus slightly distorts the CMB spectrum 
away from its Planckian shape. Seen against the CMB, clusters 
cast shadows below 217 GHz and shine at frequencies above. 
Massive galaxy clusters can be identified by this characteristic 
spectral appearance at centimetre to millimetre wavelengths if 
the angular resolution of the telescope is of order 1' or better, for 
example with the South Pole Telescope (R uhl et al.|20 04). 

The thermal SZ effect is quantified by the Compton-v param- 
eter 



m e c^ 



I 



y(s)= ~(Tt dzT(s,z)p e (s,z) , 



(1) 



which compares the mean thermal energy k%T with the rest 
energy m e c 2 of the electron and multiplies their ratio with the 
Thomson cross section crj and the electron number density p e . 
By Compton-upscattering, the specific intensity of the CMB 
seen through a galaxy cluster changes by 



A/sz(s) 



= g(x)y(s) , 



B^T) 
relative to the Planck spectrum B^T) of the CMB, where 

tiCO 

X "k^f 
is the photon energy in units of the mean thermal energy, and 



(2) 



(3) 



g(x) 



xe 



e x + l 



(4) 



describes the frequency dependence of the thermal SZ effect. 

Due to the finite telescope resolution, the observable quan- 
tity is not quite the Compton-y parameter or the specific CMB 
intensity change from Eq. Q, but rather the beam-convolved in- 
tensity change or Compton-y parameter. The Compton-y profile 
convolved with a beam profile b{s) is 



m 



2.2. Basic relations 



I 



d s'y(s')b(s - s') 



(5) 



An ideal gas assumed to have a polytropic stratification with in- 
dex y attains the density 



P = Po<P 



1/(7-1) 



(6) 



in hydrostatic equilibrium with the suitably scaled, dimension- 
less gravitational potential 



r_L^ (0cut _ 0) 



(7) 



Here, <I> is the Newtonian gravitational potential, O cut is the po- 
tential value where the density of the bound gas is supposed to 
drop to zero, and po and Pq are fiducial values for the gas den- 
sity and its pressure at an arbitrary, suitably chosen radius. By 
the ideal gas equation, the gas temperature is 



T = T <p . 



(8) 



also relative to the temperature 7o at a fiducial radius, for which 
we choose the virial radius without loss of generality. Equations 
(|6j and |9| were derived in detail in KMMSB. 

Combining Equations d6l), ((SI and (fTh, the Compton-y pa- 
rameter can be rewritten in terms of the gravitational potential <p 
as 

y(s) = ^lr po fdz^(s,z), (9) 

m s c l J 

where the exponent rj - y{y - 1) . For realistic polytropic in- 
dices, 1.1 < y < 1.2 (Finoguenov et al. 2001 ), the exponent rj is 
quite a large number, 6 < 77 < 11. 

Using Eq. p}, the beam-convolved Compton-y parameter 
can be re-written as 

y(s)= fd 2 s'y(s')b(s-s') 

" liCrT -Topofdzfd 2 s' i(P{s',z)b(s - s') (10) 
2 -T Q po I dztpi(s,z) , 



k B cr Tr 



m e c^ 



where the overbar abbreviates the convolution. 

Combining the last expression with Eq. pi leads to a result 
which can be identified with an effective, three-dimensional, pro- 
jected pressure P(r), 



g(x) jTqpq dztp 



BJT) 



n (s,z) 



I 



= dzP(s,z) 



(11) 



We shall argue now that the beam convolution can in fact 
be ignored for our purposes. Current and future SZ observations 
reach an angular resolution very much better than the angular 
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resolution that can be achieved with potential reconstructions 
based on gravitational lensing. Assuming 50 background galax- 
ies per square arc minute and averaging over 20 galaxies to ob- 
tain a sufficiently robust weak-lensing signal, the resolution of 
a weak-lensing map corresponds to » 35". Beam profiles of 
modern thermal-SZ observations are much narrower than that, 
allowing us to approximate the beam b in ( fTT| ) by a Dirac delta 
distribution. This paper uses a configuration of ALMA which re- 
sults in a beam size of 1" and is thus larger than the resolution 
of strong lensing observations (e.g. HST/ACS with a 0.1"). 
From Eq. (jTTb, we can then infer the relation 



P(r) = P <p\s,z) 



(12) 



between the effective pressure P(r) and the scaled gravitational 
potential ip, where the amplitude 



kftCTj 

Pq =g(x) yT po , 



(13) 



was introduced. 



Equation ( 12 1 serves as the basis for a method analogous to 



that presented in KMMSB for recovering the projected gravita- 
tional potential: 



1. 



By deprojection of the measured, relative specific inten- 
sity change AI^B^iT), we find an estimate for the three- 
dimensional effective pressure P. As described in KMMSB, 
this can be achieved by means of the Richardson-Lucy algo- 
rithm. 



2. Next, we use Eq. (12i to find an estimate ip for the scaled 



gravitational potential tp, 



<P = 



P(r) 
Po 



i/i 



(14) 



The estimate ip of the three-dimensional potential is then pro- 
jected to find an estimate tfr for the two-dimensional poten- 
tial, 



<Ks) 



■/ 



dzip(s,z) 



Since r\ is large, the exponent 1 jr\ is a small number, which 
is a most welcome property of Eq. ( 12 1. Fluctuations in the esti- 



mate P of the deprojected effective pressure will be substantially 
smoothed that way. 

2.3. Deprojection 

As described in detail in KMMSB, we adopt the Richardson- 
Lucy algorithm ( [Lucy 1974 1994[ > for deprojection. For sim- 
plicity, we assume spherical symmetry for now, which is not a 
necessary requirement and will be relaxed in further work. 

Generally, Richardson-Lucy deprojection connects two 
functions, say f(r) and g(s), related by the projection 



g(s) = jdrK(s\r)f(r) , 
/(r) = jdsK'(r\s)g(s), 



(16) 

(17) 



mediated by the projection kernel K(s\r) and the deprojection 
kernel K'(r\s). In our application, f(r) represents a function 
given in three dimensions and g(s) its projection onto the sky. 



For a spherically-symmetric body, the normalised projection 
kernel K(s\r) is 



K(s\r) = 



2r ©(r 2 - s 2 ) 

n y r 2 _ s 2 



(18) 



where 0(x) is the Heaviside step function. The projection and 
deprojection kernels are related by Bayes' theorem, 



f(r) 
g(s) 



(19) 



Since f(r) is unknown, so is the deprojection kernel K'(r\s). 
However, given an estimate fj(r) for the function f(r), an esti- 
mate gi(s) of the projection g(s) is 



Us) = JdrK(s\r)Mr) 
implying the estimate 



K'(r\s) = ^K(s\r) 
gi(s) 



(20) 



(21) 



for the deprojection kernel. This suggests an iterative scheme: 
Given an estimate fi{r) for the function f(r) at the iteration level 
i, an improved estimate ,//+i(r) is found by 



f M (r) 






K(s\r) . 



(22) 



Including a regularisation term for suppressing small-scale 
fluctuations, the change in the estimate fi(r) by the /-th iteration 
step is 

Af i =A H f i + A s f i , (23) 



(15) where A H fj(r) represents the change f i+i - f t given by Eq. (22 1 



and As f, may contain an entropic term S against a suitably cho- 
sen priory. 

Specialising the Richardson-Lucy algorithm to our spe- 
cific case, we replace g{s) by the observable intensity change 
AI sz B a l (T) relative to the Planck spectrum of the CMB and f(r) 
by the effective pressure P{r). 

Identifying the projected function g(s) with the radial profile 
of specific intensity change A/sze(s) and the deprojected func- 
tion f(r) with the effective pressure P{r), the complete iteration 
including an entropic regularisation term S reads 



AP, = P, 
with 



/ 



as — K(s\r) - 1 - a In + S 



AW*) 



X 



-I 



AI S z,i(s) = drK(s\ r )Pi( r ) . 



(24) 
(25) 



The complete algorithm is characterised by two parameters, 
the amplitude of the regularisation term a and the smoothing- 
scale length L, contained in the prior x- The value a — and 
L - Oh Mpc correspond to no regularisation and no smooth- 
ing, respectively. 
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Fig. 1. Azimuthally averaged and normalized specific inten- 
sity change profile of a simulated galaxy cluster with mass 
5 x 10 14 hr l M & at redshift 0.2. The displayed maximal range 
of radial values corresponds to the virial radius of the cluster. 



3. Simulating SZ observations 

3. 1 . Simulating SZ observations 

For testing the algorithm sketched above, we simulate thermal 
SZ signal of a massive galaxy cluster, assuming spherical sym- 
metry, hydrostatic equilibrium and an NFW density profile p(r) 



P(r) = 



Ps 



(r/r,)(l + r/r s y 



(26) 



with the scale radius r s - r2oo/c and the characteristic density 
p s of the halo. We choose the concentration parameter c — 5 
and a spatially flat standard ACDM cosmology with Q. m Q = 0.3, 
flbo - 0.04 and Q A o = 0.7. The properties of the ICM are chosen 
as follows: 

- The plasma contains 75% hydrogen and 25% helium by 
mass. Both components are completely ionised. 

- The gas-mass fraction equals the universal baryon mass frac- 
tion / b = Q b /n m . 

- The gas has a constant polytropic index of y — 1 .2. 

The gas density and temperature profiles are then calculated 
using Eqs. d6| and d8). To obtain a temperature profile which 
drops to zero at a large radius, we choose a large cut-off radius 
for the gravitational potential of r cut = lOO^oo- Given the den- 
sity and temperature profiles, the specific intensity change A/sz 
can be calculated using Eq. ( p"2| ). 

The scale radius of our simulated cluster is r s = 
0.25 /r'Mpc, its virial radius is r wu - 1.25 /^'Mpc. 

3.2. Background fluctuations due to unresolved clusters 

Galaxy clusters unresolved by the telescope beam contribute a 
background noise level y\, g that needs to be taken into account in 
all following calculations. 

Clusters are unresolved if they appear (much) smaller than 
the beam size. The background signal is thus dominated by 
low-mass clusters (Bartelmann 2001). Since an ideally homo- 
geneous background could be removed from the data, we only 



have to consider the average background fluctuation level Aybg- 
The mean background level contributed by unresolved clusters 
is 



Vbj 



dV 
"dl 



-/ 



(l+zf fdMn(M,z)Y(M,z) 
d 2 N(M,z) 



dVY(M,z)- 



dMdV 



(27) 



where the mass- and redshift integrations have to be carried out 
over that area in the mass-redshift plane where clusters are unre- 
solved. The integrated Compton-y parameter Y(M, z) of a cluster 
with mass M at redshift z is 



/ 



ds y(s) 



m e c 2 D 1 , 



N e 



(28) 



where the angular-diameter distance D^ to the cluster appears. 
The total number of (hot) electrons is N e . 

We choose the cluster mass function n(M, z) described by the 
Sheth-Tormen model ( |Shefh & Tormen|1999) , 



n(M,z)dAf = A 



1 + 



1 \ p dv 



• 2 « / M dM 



exp 



M , (29) 



where v = -\/aS c [(To(M)D + (z)] is the linear amplitude required 
for the collapse of a density fluctuation with the present rms fluc- 
tuation ctq(M). The linear growth factor of the density perturba- 
tions is D+ and the critical linear density contrast for non-linear 
collapse is 6 C . p is the mean background density at the present 
epoch. For the remaining parameters, She th et aT] ( |2001| l find 
A = 0.322, a = 0.707 and q = 0.3. 

The number of background sources is a discrete probability 
distribution and therefore follows a Poissonian distribution. But 
due to the very high number of background galaxy clusters the 
Poissonian distribution, according to the central limit theorem, 
converges pointwise towards the normal distribution. Therefore 
background thermal-SZ fluctuations can be assumed to be given 
by a Gaussian fluctuation of the number of clusters per unit mass 
and volume, if one also neglects any cluster correlations. Thus, 
the rms background fluctuation is 



Ay bg 



! dM ! 



, d 2 N(M,z) 

dVY 2 (M,z) ! — — 

V ' dMdV 



1/2 



(30) 



Following Eq. p0| , each point in the cluster coordinate frame 
follows a Gaussian distribution with a mean equal to the sum 
of the analytic value of the Compton->' parameter and the sig- 
nal of the unresolved clusters. The variance of the distribution is 
then equal to the rms of the background fluctuations according 
toEq. (f30i. 



3.3. Instrument noise 

Another relevant source of noise is the measurement noise of 
the telescope and the detector. We assume that the error has a 
Gaussian distribution with standard deviation cr around zero, 
with cr depending on the telescope and detector configuration 
and the setup of the specific observation. 

Since we need a resolution for our simulation that is compa- 
rable to other observations of massive galaxy clusters, we choose 
the Atacama Large Millimeter/submillimeter Array (ALMA) as 
an example telescope for achieving sufficiently precise thermal- 
SZ observations. For all simulations presented in the remainder 
of this paper, we use the following configuration of ALMA: 
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All N -32 antennae are assumed to be available with a base- 
line of 650 m. This configuration has an angular resolution of 
1.0" in a frequency band chosen to be centred on v = 116 GHz. 
Its bandwidth is assumed to be Av = 7.5 GHz. 

According to the ALMA user manual (Lundgre n|2012[ l the 
point-source sensitivity cr of ALMA in units of Jansky is 



2£b T S y S 



j] q T] c A eS y/N(N - l)«pAvf ; , 



(31) 



where T sys is the temperature of the system, j] q the quantum 
efficiency (e.g. rj q = 0.96), rj c the correlator efficiency (e.g. 
rj c = 0.88), A e ff the effective area of the antenna, n p the number 
of polarisation states (e.g. « p = 2) and ?j nt is the integration time. 
Given a required sensitivity, we choose the integration time to 
be 5 days. Given the standard deviation of the instrument noise, 
the final observed image can be simulated assuming a Gaussian 
distribution of the signal with the variance a 2 . 

For this configuration, the radial profile of the thermal SZ 
signal can be obtained by averaging the simulated image in cir- 
cles around the cluster centre. This profile is shown in Fig. [Tlfor 
one realisation of the simulated galaxy cluster. At a radius near 
0.8 /!~'Mpc (i.e. 3.5r s ), the cluster signal sinks below the instru- 
mental noise and the background-fluctuation level. 




0.4 0.6 0.6 

r/Mpc/h 



Fig. 2. Reconstructed scaled gravitational potential tp of a simu- 
lated galaxy cluster with mass 5x 10 14 h~ l M Q at redshift 0.2. The 
potential was reconstructed using a = 0.2 and L - 0.3 h Mpc 
in the regularisation term. 



4. Results 

Applying our algorithm to the simulated specific intensity 
change A/sz returns an estimate for the three-dimensional effec- 
tive pressure P. Through Eq. (J6l), an estimate of the gravitational 
potential <p can then be obtained. For the specific data set shown 
in Fig.^| the estimate for the gravitational potential is shown 
in Fig. p[ reconstructed with a — 0.2 and L - 0.3 h~ x Mpc. As 
expected, the reconstruction is noisy due to the low signal-to- 
noise ratio in the input data themselves, even though taking the 
azimuthal average significantly increases the signal-to-noise ra- 
tio. The reconstruction remains reliable to a maximum radius of 
r^O.S^MpcCcf.Fig.Q, 



4.1. Testing the algorithm 

As we have shown in KMMSB, the Richardson-Lucy deprojec- 
tion algorithm works reliably for reasonable amplitudes a of the 
regularisation and smoothing scales L. We also showed there that 
an increasing regularisation amplitude can substantially suppress 
the fluctuations in the reconstructed X-ray surface brightness 
(see Fig. |3t in KMMSB). Plots showing analogous results for 
the potential reconstruction from thermal SZ data are shown in 
Fig.pa for < a < 0.9 with a constant L - 0.3 hr l Mpc. Varying 
the smoothing scale while keeping a constant leads to qualita- 
tively similar results. For large amplitudes a, the fluctuations 
in the recovered potential are damped considerably. For large 
smoothing scales L, the noise is overestimated at large radii. 

The effect of varying a and L on the recovered lensing po- 
tential is marginal for radii below r « 0.8 /r'Mpc. Only at radii 
larger than that, the lensing potential becomes smoother and 
fluctuates less, but tends to overestimate the true lensing poten- 
tial due to the normalisation conditions of the algorithm. Again, 
changing L while keeping a fixed has qualitatively similar ef- 
fects. 




Fig. 4. Relative rms deviation of the projected potential iff from 
its true profile i^ true , derived from 200 realisations of the mod- 
elled galaxy cluster with mass 5 x IQ U h~ l M Q at redshift 0.2. 
The blue line shows the 5 % level for comparison. 



4.2. Error estimation 

We quantify the error of our algorithm by 200 bootstrap resam- 
plings of the galaxy cluster model described above, to which we 
applied our reconstruction algorithm. From all reconstructions, 
the rms deviation of the recovered projected potential from its 
true profile is then calculated, 



rms(s; iff, i/w) = 



1 N 

-Y 



(^^S) - Cue m (*)) 



C°ue m (*r 



1/2 



(32) 



where quantities with a superscript 'norm' are normalised to 
reach zero at the maximum projected radius. 

The result of this bootstrap is shown in Fig. HI together with 
a reference line at the 5 % level of deviation (blue dotted line). 
We achieve a relative accuracy of less than 5 % for radii smaller 
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Fig. 3. (a) Comparison between potentials recovered with different choices for the amplitude a of the regularisation function in 
the reconstruction of the effective pressure P, keeping the smoothing scale constant at L - 0.3 h Mpc. The result expected from 
Eq. (|6| is plotted for reference. Varying the smoothing scale within 0.1 < L < 0.9 h~ l Mpc leads to qualitatively similar results, (b) 
As (a), but for the projected potential i//. 



than as 0.9 /r'Mpc. At larger radii, the rms reaches values of up 
to < 15 %. This is clearly due to the increasing noise in the signal 
at large radii. 

5. Conclusions 

Following the X-ray analysis presented in the first paper of this 
series, we have shown in this paper how the observable pro- 
vided by the thermal SZ effect in clusters, i.e. the relative inten- 
sity change of the CMB observed through the hot intracluster 
plasma, can be converted into the projected, two-dimensional 
gravitational cluster potential. As in the preceding paper, the 
goal of this study is to bring all cluster observables - strong 
and weak gravitational lensing, X-ray emission, the thermal SZ 
effect and ultimately also galaxy kinematics - on a common 
ground to use all of them in a joint reconstruction procedure 
recovering the gravitational potential best compatible with all 
these observables. 

Assuming hydrostatic equilibrium between the hot gas and 
the gravitational potential, and further assuming a polytropic gas 
stratification, we have derived how the Compton-y parameter re- 
lates to the gravitational potential. This allowed us to construct 
an algorithm beginning with the Richardson-Lucy deprojection 
of the observed, two-dimensional thermal-SZ intensity change 
into the three-dimensional, effective pressure proportional to the 
Compton-y parameter. Richardson-Lucy deprojection is the first 
step in the algorithm requiring symmetry assumptions. For sim- 
plicity, not by necessity, we have chosen to assume spherical 
symmetry for this initial study. The deprojected Compton-y pa- 
rameter is then readily converted to the three-dimensional grav- 
itational potential, which can finally be projected. 

Our implementation of the Richardson-Lucy deprojection al- 
gorithm contains an entropic regularisation term with two pa- 
rameters, an amplitude and a smoothing scale. Both suppress 
the reconstruction noise as they should, rendering the result- 
ing two-dimensional potential slightly dependent on their values. 
For quite wide ranges of reasonable parameter choices, however, 
the result is very close to the expected, two-dimensional cluster 



potential known from the cluster model underlying the simula- 
tions. 

We have tested this algorithm with synthetic thermal-SZ 
data simulated with a spherically-symmetric body of mass 5 x 
10 14 h~ l M at redshift 0.2, supposed to be observed with the 
signal-to-noise characteristics of one specific configuration of 
the ALMA interferometer. In addition to instrumental noise, we 
have included background fluctuations in the thermal SZ signal 
due to unresolved clusters. The results look very promising: The 
three- and two-dimensional gravitational potentials are very well 
reproduced. Bootstrapping shows relative rms accuracies of the 
recovered, two-dimensional potential at or below the 5 % level 
can be achieved at cluster-centric radii r < 0.9 h~ l Mpc. 
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